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Abstract 

The vibration and performance opti- 
mization procedure using regression ana- 
lysis has been successfully applied to an 
advanced aeroelastic blade design study. 
The major advantage of this regression 
technique is that multiple optimizations 
can be performed to evaluate the effects 
of various objective functions and con- 
straint functions. In this application, 
the data bases obtained from the rotor- 
craft flight simulation program C81 and 
Myklestad mode shape program are analy- 
tically determined as a function of each 
resign variable. Those predicted results 
from regression equations, such as per- 
formance, vibration, and modal parameters, 
when compared with C81 and Mykle 3 tad out- 
puts, correlate exceptionally well. The 
regression equations also predicted the 
minimum of 4/ rev total vertical hub shear 
based on the coefficients of each equa- 
tion. This approach has been verified for 
various blade radial ballast weight loca- 
tions and olade planforms. This method 
can also be utilized to ascertain the 
effect of a particular cost function which 
is composed of several objective functions 
with different weighting factors for 
various mission requirements without any 
additional effort. Utilization of this 
technique can significantly reduce the 
engineering efforts and computer time to 
optimally design a high performance and 
low vibration blade. 

Introduction 

It is highly desirable for most heli- 
copter engineers to design a vehicle 
having high performance and low vibrations 
(References 1-13). With a best dynamic 
blade as an input to the airloads program, 
the blade having minimum vibration and 
maximum performance under certain con- 
straints could be determined by using an 
existing optimization codes or vice versa, 
from an optimized airloads distribution to 
find a desired blade planform. Blade 
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dynamic and aerodynamic effects are 
coupled within the design range of inte- 
rest. Separation of these effects during 
the design procedure may not be possible 
to obtain the best result that one ex- 
pects. Therefore, the approach which can 
be utilized to optimize dynamic and aero- 
dynamic effects is stro: jly recommended. 

Vibration and performance data gene- 
rated from C81 and the coefficients of 
modal participation factor (CMPF) of hub 
shear and hub moment generated from Mykle- 
stad can be analytically expressed as a 
function of each design variable using re- 
gression analysis (References 14 - 20) . 
Regression equations not only directly 
provide the sensitivity of each blade de- 
sign variabla, but also combine both dy- 
namic and aerodynamic effects within the 
overall design procedure. Furthermore, 
repression technique need not be performed 
in a continuous run; it may be carried out 
individually or in groups, as convenient. 
This technique can also treat numerous 
design variables, objective functions, 
constraint functions, and various combina- 
tions of several objective functions in a 
convenient manner. After the data base is 
obtained from tne technique program, the 
optimization criteria can be varied, based 
on various mission requirements. There- 
fore, a significant savings on computer 
time and engineering efforts have been 
achieved. 

The optimization procedure of the re- 
gression analysis was first used at Kaman 
in its analytical studies of the Control- 
lable Twist Rotor in developing secondary 
control requirements to minimize vibra- 
tion, with constraints on horse-power, 
angle of attack, and blade bending mo- 
ments. This control optimization was done 
for both steady and one-per-rev controls, 
as well as for higher harmonic controls 
(Reference 20). Blade controls on the 
full scale Multicyclic Controllable Twist 
Rotor with higher harmonics were optimized 
experimentally by using wind tunnel re- 
sults for the data base (Reference 17). 

The optimization procedure was also 
used to investigate the effects of several 
blade design parameters as Independent 
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variables in a study of an advanced flight 
research rotor (Reference 18) . All pre- 
vious results are obtained either from 
Raman' s program, 6F, or from the wind tun- 
nel test, on the hinged blade. The input 
mode shapes for the 6F are uncoupled modes 
with pitch horn control and servo flap 
control degrees of fr Mom. 

The main rotor this study is a 
hingeless, 4-bladed General Purpose Re- 
search Rotor (GPRR) (References 18, 19) 
which weighs 287,5 lbs per blade and has 
27 ft radius, 25.5 in. thrust weighted 
chord, 256 rpm angular speed, and 723.8- 
ft/sec tip speed. Bingham RC airfoil 
tables are used to determine blade aero- 
dynamic coefficients. The fuselage has 
18,400 lbs total gross weight and 23 
square ft flat plate drag area. C81 was 
modified to incorporate variable sweep 
stations along the blade radial direction. 

Thirty-six C81 quasi-static (QS) trim 
cases as a function of blade built-in 
twist, sweep angle, percent tip taper, and 
taper ratio have been generated to find 
the regression equations for prrformance 
analysis at five different airspeeds from 
hover to 160 knots. The predictions of 
the horsepower from regression equations, 
which are not included in those 36 QS trim 
cases, compared with C81 are within 1.5Z 
of the total range of interest. 

The regression equations of the modal 
parameters also have been generated using 
84 Myklestad cases by adding blade ballast 
running weights along the blade radial 
stations. The predicted results from re- 
gression equations, compared with the 
Myklestad, are in excellent agreement up 
to the first six modes. 

Thirty-five C81 QS Trim, followed by 
Time-Varying Trim (TVT) , cases as a func- 
tion of blade built-in twist, percent tip 
taper, and taper ratio are used for vibra- 
tion analysis, The multiple correlation 
factors for horsepower, 4/rev vertical hub 
shear, oscillatory beamwise and chordwise 
bending moments, and torsional moments are 
correlated at least 95. 4Z. The excellent 
predictions from regression equations for 
the vibration data are also presented. 

With the exceptionally well-fitted 
regression equations from C81 and Mykle- 
stad, the blade can be dynamically con- 
trolled by controlling each individual 
CMPF , or its product with MPF, to achieve 
the design goal under certain constraints. 

Performance Analysis 

In order to determine blade charac- 
teristics for the performance analysis, 
blade physical parameters from Reference 


18 are used for the baseline blade. Blade 
torsional control spring is pre-determined 
as an input to Myklestad coupled mode 
shapes program such that the blade clamped 
torsional frequency is 25.6 Hz. Only the 
first out-of-plane mode shape is used as 
an input in C81 for the performance study, 
Bingham Rfc 8Z, 10Z, and 12Z airfoil tables 
are used to look up blade local aerody- 
namic lift, drag, and pitching moment co- 
efficient's. Blade built-in twist, sweep 
angle, percent tip taper, and taper ratio 
are treated as independent variables input 
to C81 to vary blade airloads distribu- 
tion. QS trim in C81 uses the first flap- 
ping mode; therefore, blade sweep angle 
gives no dynamic coupling effects t.nd only 
has aerodynamic effects on Mach No, reduc- 
tion; and aerodynamic effects on pitching 
moment variation due to aerodynamic center 
shift. The blade sweep station starts at 
that point at which the Mach No. is the 
same as the Mach No. on the blade tip. 
There are four independent variables in 
the analysis, and the range of interest of 
these variables is listed in Table 1. 


Table 1. Independent variables for 
performance analysis. 


Independent 


Variables 

Levels 

Built-in Twist* 

-8°, -12°, -14®, -16° 

Sweep Angle** 

20°, 0°, -20® , -30® 

Percent Tip Taper 

15:, 25Z, so: 

Taper Ratio 

1.1:1, 2:1, 3:1 

*Built-in Twist; + Nose Up 

[ **Sweep Angle; 

+ Forward Sweep 


The quadratic regression equation of the 
independent variables is written as 
follows : 


A ° + Ji Vi * ih a ‘i‘i 


N-l 


+ l I 

i-1 j-i+1 


A ij 6 i S j 


Where Y is the dependent variable; 6 is 
the independent; and A is the coefficient 
of regression equation. 


There are 144 different combinations 
for these variables. Only 36 combinations 
are randomly selected as inputs for C81 QS 
trim at each flight speed. The regression 
equations having linear and quadratic 
terms which are generated from these data 
are shown in Table 2 for 5 different 
speeds. These equations give multiple 
correlation coefficients of 97. 5Z or 
better at each different speed from hover 
to 160 knots. With the existing data 
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Table 2. Regression equations for performance analysis at 5 different flight speeds. 


Coefficient 

Variable 

Hover 

40 Knots 

80 Knots 

120 Knots 

160 Knots 

A o 


1975.08 

1324.90 

1038.82 

1316.16 

2054.82 

*1 

4 i 

" 

0.85 (9) 


-0.45 (7) 

-1.31 (7) 

a 2 

f 2 

12.13 (5)* 

17.96 (4) 

18.89 (7) 

32.23 (3) 

44.88 (3) 

A 3 

4 3 

-77.59 (7) 

-47.36 (8) 

-27.60 (8) 

-28.33 (10) 

-27.88 (4) 

A 4 

4 4 


— 

-11.31 (3) 

20.65 (5) 

- “ - 

A 11 

V«1 

-0.067 (2) 

-0.04 (2) 

-0.04 (1) 

-0.05 (1) 

-0.08 (1) 

*22 

V 6 2 


0.52 (6) 

0.74 (6) 

1.46 (2) 

2.078 (2) 

A 33 

*3* 6 3 

19.43 (6) 

10.91 (7) 

5.45 (9) 

2.84 (11) 

- - - 

A 44 

V 4 4 


130.55 (11) 

63.80 (10) 

106.75 (12) 

140.25 (9) 

A] 2 

4 1* 4 2 

- - - 

0.04 (10) 

-0.01 (4) 

-0.05 (6) 

-0.10 (6) 

A 13 

4 1* 4 3 

— 

-0.11 (12) 


-0.08 (13) 

- - - 

A 14 

4 1* 4 4 



... 

0.43 (14) 

0.61 (10) 

A 23 

4 2* 4 3 

1.26 (1) 

0.10 (1) 

-0.29 (5) 

-1.39 (9) 

-1.84 (5) 

A 24 

V 4 4 

-7.0 (3) 

0.75 (5) 


6.68 (8) 

9.12 (8) 

A 34 

4 3* 4 4 

-106.88 (4) 

-77.88 (3) 

-38.14 (2) 

-22.17 (4) 


M.C.C.** 


0.987 

0.991 

0.982 

0.975 

0.978 

S.E.E.*** 


11.7 

5.8 

4.7 

7.2 

9.1 

6. Sweep 

Built-In Twist 6, 

. Taper Ratio 

6^ % Tip Taper 


1 * Sensitivity ** Multiple Correlation Coefficient *** 

Standard Error of the Estimate 



Table 3. Regression equation for performance analysis with airspeed as an 
independent variable. 


Coefficient 

Variable 

Horsepower 

" ' 

Coefficient 

Variable 

_ _j 

Horsepower 

Aq 

A 1 

4 i 

2063.68 

A 12 

4 l‘ 4 2 


A 2 

4 2 

36.55 (8)* 

A 13 

4 1* 4 3 


A 3 

4 3 

-21.58 (4) 

A 14 

4 1* 4 4 


A 4 

4 4 

- - - 

A 15 

4 1* 4 5 


A 5 

4 5 

-22.76 (2) 

A 23 

4 2* 4 3 


A 11 

V 4 i 

-0.06 (3) 

A 24- 

4 2* 4 4 

2.01 (6) 

A 22 

Y 4 2 

1.06 (11) 

A 25 

4 2* 4 5 

-0.12 (7) 

A 33 

4 3* 4 3 

- - - 

A 34 

4 3* 4 4 

-57.54 (10) 

A 44 

V«4 

- - - 

A 25 

4 3* 4 5 

0.23 (5) 

A 55 

V«5 

0.13 (1) 

A 45 

4 4* 4 5 

1.03 (9) 


Multiple Correlation Coefficient: 

0.999 


Standard Error of the Estimate: 

18.4 


«j Sweep 4 3 

Taper Ratio 

Airspeed 

<2 Built-In Twist 4^ 

% Tip Taper 

* Sensitivity 
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base, the regression equation for perfor- 
mance as a function of aiiapeed has also 
been analyzed. The multiple correlation 
coefficient from the equation with air- 
speed as an independent variable is cor- 
related at 99 . 91 , shown in Table 3. 

The sensitivity results from re- 
gression equations show that each design 
variable has a clear performance trend at 
each airspeed and for the airspeed sweep. 
The independent variables in these regres- 
sion equations have not been normalized. 
Therefore, the physical parameters are 
treated as the input to these regression 
equations. From Table 2, blade sweep 
angle squared, built-in twist squared, and 
built-in twist are the three most impor- 
tant terms at 160 knots from the per- 
formance regression equation sensitivity 
result. Also, the product of taper ratio 
and built-in twist, sweep angle squared, 
the product of built-in twist and percent 
tip taper, and the product of percent tip 
taper and taper ratio are the four most 
important terms in hover. Blade sweep 
angle squared, the product of taper ratio 
and percent tip taper, percent tip taper 
and the product of sweep angle and built- 
in twist are the four most important terms 
in the regression equation at 80 knots. 
From Table 3, the regression equation 
shows that airspeed squared, airspeed, 
blade sweep angle squared, and taper ratio 
are the four most important terms in the 
whole airspeed sweep region. Also from 
Table 2, the regression equation shows 
that the constant term has the minimum 
value at 80 knots. All the design vari- 
ables have either an increased or a de- 
creased contribution to the constant term 
at each flight speed, depending on the 
combination of each individual design 
variable. 

In order to gain a better under- 
standing of the effects each independent 
variable contribution to performance, the 
plots of horsepower vs each Independent 
variable at different speeds (Fig. 1 to 4) 
are described as follows: 

1. For a blade having -10° built- 
in twist and 25Z tip taper, results show 
that a 3:1 taper ratio blade saves 20 HP 
over a 1.1:1 taper ratio blade at 160 
knots i saves 25 HP at 80 knots; and saves 
80 HP in hover. 

2. Results also show that a 30° aft 
sweep blade saves 75 HP and 35 HP at 160 
knots, 60 HP and 35 HP in hover, and 35 HP 
and 25 HP at 80 knots over a non-swept 
blade and a 20° forward sweep blade, 
respectively. 

3. The 3:1 tapet ratio blade has 
better performance than the 2:1 and 1.1:1 


taper ratio blades at all flight speeds of 
interest, except at 160 knots. 

4. For the blade having -16° built- 
in twist and 50Z tip taper, the 3:1 taper 
ratio blade uses slightly more HP than a 
1.1:1 taper ratio blade at 160 knots, and 
saves 150 HP in hover and 40 HP at 80 
knots . 

5. For a high negative built-in 
twist blade, -16° , the best performance 
is at hover, with very little effect on 
performance at 160 knots. The best per- 
formance at 160 knots is with the blade 
which has approximately -10° built-in 
twist. 

The prediction of the horsepower from 
regression equations compared with C83 
trim results is exceptionally good. The 
difference between the two result? is 
within 1.5Z of the total range of inter- 
est. The comparison is shown on Tables 4 
and 5. 


The regression equations for horse- 
power at 160 knots, 80 knots, and hover 
are used for the performance optimization 
study. Power limits from C81 QS trim are 
treated as constraints at 160 knots and 80 
knots. Those constraints for maximum 
power available are assumed to be 1740 HP 
at 160 knots and 840 HP at 80 knots. The 
minimum horsepower from 36 QS trim cases 
used as the starting point for optimi- 
zation is the blade having a plan, rrn 30° 
aft sweep, -16° built-in twist, 3:1 taper 
ratio, and 50Z tip taper. The optimiza- 
tion code KAOPT (Reference 21) is used for 
performance optimization. There are two 
minimum points detected using the KAOPT 
volume search technique. The first point 
is the blade having $ 0° aft sweep, -15.8° 
built-in twist, 50Z tip taper, and 3:1 
taper ratio. The second point is 20° for- 
ward sweep, -10.4° built-in twist, 44Z 
tip taper, and 3:1 taper ratio. The per- 
formance results are 1740 HP, 822 HP, 
1500 HP for point 1; and 1740 HP, 841 HP, 
and 1616 HP for point 2 at 160 knots, 
80 knots, and hover, respectively (also 
shown in Table 4). The contour plots of 
power at hover, with and without con- 
straints, are shown in Fig. 5. For 
1740 HP available constraint applied to 1 
g thrust, 160 knots and 1.5 g thrust, 120 
knots, level flight conditions, the mini- 
mum power at hover within constraints is 
1516 HP, and the blade has 30° aft sweep 
-14.54° built-in twist, 3:1 taper ratio, 
and 50Z tip taper. 


Modal Analysis 


The elastic rotor usee seven indepen- 
dent modes representation in the C81 air- 
loads analysis. The time history of rotor 
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Table 4. Performance predictions for regression equation vs C81 at 
three different speeds. 


4 1 

*2 


4 < 

V » 160 

KNOTS 

V = 80 KNOTS 


V - 0 KNOTS 

SWEEP 

BUILT 

IN 

TAPER 

l TIP 

REGRESSION 

C81 

REGRESSION 

C81 

REGRESSION 

C81 

TWIST 

RATIO 

TAPER 

(HP) 

(HP) 

(HP) 

(HP; 

(HP) 

(HP) 

- 20* 

- 13° 

2.5:1 

301 

1753.65 

1776.25 

849.69 

1 ^ 

1624.21 

1618.44 

- TO* 

- 9* 

2.5:1 

301 

1712.66 

1715.82 

838.76 

839.81 

•643.57 

1629.31 

- n* 

- 9.4" 

3:1 

461 

1766.20 

1764.91 

851.47 

847.61 

1624.22 

1627.86 

- 30° 

- 11.8° 

3:1 

461 

1698.0 

1709.82 

816.0 

814.16 

1559.0 

1540.03 

- 30° 

- 13* 

2.5:1 

301 

1/11.29 

1714.56 

830.20 

829.73 

1590.81 

1561.65 

- 20* 

- 9* 

2.5:1 

301 

1750.94 

1752.02 

857.88 

861.98 

1676.96 

1678.20 

20° 

- 13" 

2.5:1 

301 

1761.53 

1768.76 

864.52 

852.3? 

1624.21 

1618.01 

20° 

- 9* 

2.5:1 

301 

1742.53 

1748.47 

861.22 

862.08 

1676.96 

1677.67 

- 30* 

- 15.8* 

3:1 

501 

1740.0 

1743.76 

822.0 

020.66 

1500.0 

1508.61 

- 30* 

- 13" 

1.5:1 

301 

1715.24 

1716.92 

843.64 

839.34 

1639.15 

1616.54 

- 30* 

- 9" 

1.5:1 

301 

1723.97 

1724.64 

853.37 

848.83 

1686.86 

1662.44 

- 20" 

- 13" 

1.5:1 

301 

1757.60 

1770.04 

663.13 

865.69 

1672.55 

1662.93 

- 20° 

- 9* 

1.5:1 

301 

1762.26 

1760.44 

872.50 

875.68 

1720.26 

17' 1 0.64 

20* 

- 10.4" 

3:1 

441 

1740.0 

1749.87 

841.0 

839.35 

1616.0 

1604.04 

20* 

. r 

1.5:1 

301 

1753.85 

1755.93 

875.84 

875.93 

1720.26 

1719.92 

20* 

- 13" 

1.5:1 

301 

1765.49 

1767.46 

867.96 

865.14 

1672.55 

1662.34 


Table 5. Performance predictions for regression equation vs C81 with airspeed 
an Independent variable. 


4 i 

4 ? 

*3 

«« 

4 S 

HORSEPOWER 


BUILT 

IN 

TAPER 

% TIP 

;*rspeeo 

REGRESSION 

C81 

SWEEP 

TWIST 

RATIO 

TAPER 

(KNOTS) 

(HP) 

(HP) 

- 30* 

• 9.0* 

1*5:1 

30% 

160.0 

1716.38 

1724.64 

- 20" 

• 13.0* 

2.5:1 

30% 

160.0 

1765.17 

1776.25 

- 16* 

- 13.6* 

3:1 

44% 

147.0 

1507.30 

1518.51 

12* 

- 14.4* 

3:1 

467 

80.0 

845.81 

850.0? 

- 11* 

- 9.4* 

3:1 

46% 

114.0 

1034.92 

1023.69 

20* 

- 10.4* 

3:1 

*41 

0.0 

1624.73 

1604.04 

- 20* 

- 9.0* 

1.5:1 

30% 

0.0 

1733.29 

1720.64 

- 30* 

- 15.8* j 

3:1 

50% 

160.0 

1761.64 

1743.76 

12* 

• 14. 4^ 

3:1 

46% 

57.0 

887.20 

880.20 

- 20* 

- 13.0* 

1.5:1 

30% 

0.0 

1677.99 

1662.93 

- !•• 

- 13.6* 

3.: 

44% 

111.0 

1001.60 

1013.50 

20* 

- 9.0* 

2.5:1 

30% 

80.0 

875.04 

862.08 

• u* 

- 13.6* 

3:1 

44% 

63.0 

8C9.24 

853.93 

20* 

- 13.0* 

1.5:1 

30% 

160.0 

1766.93 

1767.46 

12* 

- 14.4* 

3:1 

46% 

137.0 

1344.81 

1357.36 

- 30* 

* 11.8* 

3:1 

46% 

160.0 

1719.79 

1709.62 

• 11* 

• 9.4* 

3:1 

46% 

97.0 

909.68 

900.60 
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Fig. 1. Regression equation - performance vs sweep angle at 160 knot 
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r ig. 2. ^egression equation - performance vs sweep angle at 80 knots, 
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Fig. 3. Regression equation - performance vs sweep angle at hover. 
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Fig. 4. Regression equation - performance vs built-in twist at 160 knots. 
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hub shear and hub tmraent at any given sta- 
tion can be computed from the modal parti- 
cipation factors (MPF) for the last rotor 
revolution in C81. Multiply the MPF for 
each given mode by the hub shear or the 
hub moment coefficient of that mode at any 
station and sum over all modes to get the 
value at that time point. Those coeffi- 
cients of MPF can be obtained from the 
Myklestad coupled mode shape program. 
Regression analysis can be used to tune 
the coefficient of the modal participation 
factor (CMPF) or its product with MPF for 
aeroelastic blade design technique. 

The baseline blade is divided into 
nineteen 13-inch-lung equal segments, with 
each segment treated as an independent 
variable in the regression modal analysis. 
The regression equations of the first 
three out-of-plana (OP) frequencies, 
second and third OP deflections, static 
moment, flapping inextia, Lock number, and 
the CMPF of hub shears and moments of the 
first seven independent modes have been 
generated by adding blade ballast running 
weights of 1, 2, or 3 lb/in., with a total 
constant ballast weight of 39 lbs on each 
baseline blade. 

There are 6,859 possible combinations 
for putting ballast weight in a blcde with 
19 independent variables and 3 1 ' els for 
each. 84 cases are randomly o xected to 
provide enough data for linear regression 
analysis. The linear regression equation 
with 19 independent variables is written 
as follows: 

19 

Y - A q + l A i4i 
i-1 

The out-of-plane components of the CMPF of 
hub shear and moment have been curve fit- 
ted up to 7 independent modes based on a 
1-inch tip deflection, or 10° tip torsion. 
Since CMPF of hub shear of the first in- 
plane mode is either 0 or 1 , from Mykle- 
stad, no regression analysis is needed for 
that mode. At least 250 more cases are 
required if quadratic regression equations 
are considered in the modal analysis. 

The regression equations for the 
modal analysis are shown in Table 6. The 
multiple correlation coefficients (MCC) 
from the regression equations are ext- 
remely well-fitted and correlated from 
94, jZ to 99. 9Z for Myklestad modal data. 
For the first three OP frequencies, MCC 
correlates those frequencies from Mykle- 
stad output at least 98. 7Z. For static 
moment, flapping inertia, and Lock no., 
the MCC correlates no less than 99.7 Z. 
The mode shape deflections of second 


and third modes correlate better than 
98. 5Z. For CMPF of hub shear and moment 
of the first four OP modes, the MCC cor- 
relates at least 96. 1Z, and correlates 
torsion mode at least 95. 1Z. 

The regression equation sensitivity 
results are also concluded as follows: 

1. Blade outboard stations 16, 17, 
18, and 19 are very sensitive to the first 
three OP frequencies. The intercepts of 
these OP frequencies are 1.0896 P, 2.5074 
P, and 4.5889 P, respectively. 

2. Adding ballast weight in these 
four stations (16, 17, 18, 19) will de- 
crease the first OP frequency and increase 
the second and third OP frequencies. How- 
ever, adding weight at the first blade 
station will increase the first three OP 
frequencies; and adding weight at station 
8 will decrease first and second OP fre- 
quencies and increase the third OP 
frequency. 

3. The values of static moment and 
flapping inertia are increased by adding 
ballast weight in blade outboard stations 

18 and 19. However, by adding weight at 
inboard blade stations 1, 2, and 3, these 
values are decreased. Reverse trend is 
obtained for Lock number by adding the 
same ballast weight at the same stations. 

4. For the second and third OP mode 
shape deflections, putting ballast weight 
at stations 18 and 19 will make minimum 
deflections of these modes more negative 
and maximum deflection of the third OP 
mode more positive. However, adding bal- 
last weight at stations 11, 12, 13, and 14 
gives the reverse trend of the second and 
third OP modes minimum deflections and the 
same trend of the third OP mode maximum 
deflection. 

5. The CMPF of hub shear and moment 
of the first OP mode are decreased by 
adding blade ballast weight. Adding bal- 
last weight at stations 17, 18, and 19 
gives the second and fourth OP mode CMPF 
of hub shear and moment more negative and 
the 1st torsion mode less negative. Also, 
adding ballast weight at stations 18 and 

19 increases the CMPF of hub shear and 
moment of the third OP mode. 

The predicted results from the re- 
gression equations, compared with the 
Myklestad, are extremely well as shown on 
Table 7. The first three out-of-plane 
frequencies, static moment, flapping 
inertia, and Lock no. are within 1Z. The 
predicted second and third OP deflections 
are within 2.5Z. The predicted coef- 
ficients of hub shear and moment for the 
first 6 independent modes are in excellent 
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agreement, except for the fourth OP hub 
shear and moment. Higher order terms in 
the regression equation are required in 
order to have better prediction for the 
modal parameters higher than the seventh 
mode. However, the seventh mode, or 
higher than seventh modes, normally gives 
very little effect on blade performance 
and vibration analysis; therefore, the 
linear regression analysis is an appro- 
priate approach for future blade design 
study. 

Vibration Analysis 

Thirty-five C81 QS Trim followed by 
Time-Varying Trim (TVT) cases , each having 
two ballasting configurations, as a func- 
tion of blade built-in twist, percent tip 
taper, and taper ratio are used for vi- 
bration analysis at 160 knots. 

During the TVT, only flapping angles 
of the time-variant rotor are allowed to 
vary; the fuselage and control positions 
are held fixed at the values determined by 
the QS trim. The hub shear, hub moment, 
horsepower, and modal participation factor 
are obtained after the rotor reaches 
steady state within 8 rotor revolutions. 
Linear and quadratic terms are adapted to 
determine regression equations for horse- 
power, 4/rev vertical hub shear, blade 


root oscillatory beamwise and chordwise 
bending moments, and torsional moments. 

The coefficients of the regression 
equations and the multiple correlation 
coefficients are shown on Table 8. The 
multiple correlation coefficients for 
horsepower, 4 /rev vertical hub shear, 
blade root oscillatory beamwise and chord- 
wise bending moments, and torsional mo- 
ments are correlated at least 95. 4Z. 

The predictions between regression 
equations and C81 TVT results are shown in 
Table 9. The prediction of performance, 
bending moment and 4/rev vertical hub 
shear are correlated very well with C81 
TVT and the regression equation results. 

The best performance blade obtained 
from the regression equation prediction is 
a 3:1 taper ratio, -10° built-in twist, 
and 50Z tip taper blade. The 1.1:1 taper 
ratio blade has lower 4/rev total vertical 
hub shear than those blades which have 2:1 
or 3:1 taper ratio at 160 knots, from the 
regression analysis. 

Also, three different planforms com- 
bined with various ballasting configura- 
tions along the blade span have been 
investigated. There are twelve different 
ballast weight locations chosen from the 


Table 7. Regression equation prediction vs Myklestad for modal results. 
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Table 8. Regression equation for C81 vibration results 
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Table 9. Regression equation prediction vs C81 TVT results. 
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no ballast blade by putting the ballasting 
at maximum and minimum deflections and 
nodal points of the OP modes. 

The regression equations obtained 
from the combination of CMPF from Mykle- 
stad and MPF from C81 provides the sensi- 
tivity of each design variable, and also 
predicts two local minimum points of 4 /rev 
total vertical hub shears from the coef- 
ficients of each equation, shown in Fig. 6 
and 7. From these figures, the inboard 
minimum 4/rev vertical hub shear ballas- 
ting location is between station 129 and 
155, and the outboard minimum hub shear 
ballasting location is between station 246 
and 272. 

Because the GPRR blade has a large 
third OP modal component contribution to 
the 4/rev total vertical hub shear from 
the modal analysis, the inboard ballasting 
location does not have strong coupling 
between modal forces and mode shapes. 
Therefore, the results show tha the best 
vibration and performance blades for each 
of the three inboard ballasting config- 
urations have converged to the same blade 
planforms, respectively, for each ballas- 
ting location. For the inboard converged 
point, the unlapered blade predicts a 
higher power requirement and less vibra- 
tion, compared with tapered blades. How- 
ever, this trend is reversed for the 502 
tip taper, 3:1 taper ratio, and -10° 
built-in twist blade. 

For the outboard minimum, the data 
shows that there is a strong modal force 
and node shape coupling which signifi- 
cantly reduce the third OP modal compo- 
nents. For the outboard minimum point, 
the best performance blade has a similar 
vibration level compared with the unta- 
pered blade, but the performance is 152 
better than the untapered blade. Further 
study is required to investigate other 
possible local minimum vibration locations 

With the exceptionally well-fitted 
regression equations from C81 and Mykle- 
stad, the blade can be dynamically con- 
trolled by controlling each individual 
CMPF, or its product with MPF, to achieve 
the design goal under certain constraints. 
The best performance blade, obtained from 
the best ballasting configuration in this 
study, has at least 2.5 times the reduc- 
' ion of vibration level compared with the 
o/iginal ballasting configuration with 
various planform and the power requirement 
is at least 152 better than the untapered 
blade. 

Conclusions 

From the performance, modal, and 
vibration analysis of the advanced 


aeroelastic blade design study, the fol- 
lowing conclusions can be obtained from 
the results: 

1. With the exceptionally well- 
fitted regression equations from C81 and 
Myklestad, regression technique can be 
used for vibration analysis, modal anal- 
ysis, and performance analysis for de- 
signing future advanced aeroelastic rotor 
blades . 

2. Multiple optimizations can be 
performed to evaluate the effects of 
various objective functions and constraint 
functions, or to evaluate the combinations 
of several objective functions with dif- 
ferent weighting factors for various mis- 
sion requirements. 

3. Regression technique can di- 
rectly determine the sensitivity of each 
blade design variable and analyze the dy- 
namic and aerodynamic effects during the 
entire design process. 

4. The predicted results from re- 
gression equations for performance analy- 
sis, modal analysis, and vibration analy- 
sis are exceptionally good when compared 
with C81 and Myklestad outputs. 

5. For the GPRR blade, the combina- 
tion of CMPF from Myklestad and MPF from 
C81 predicts the same converging points 
for different blade planforms and differ- 
ent ballast weight configurations along 
the blade. 

6. The best performance blade ob- 
tained from the best ballasting configura- 
tion has at least 2.5 times the reduction 
of vibration level when compared with 
original configurations and the power 
requirement is at least 152 better than 
the untapered blade. 
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OPTIMAL DESIGN APPLICAflON ON THE ADVANCED AEROELASTIC ROTOR BLADE 

Fu-Shang Wei 
and 

Robert Jones 


Dave Peters. Washington University : Do you have some feel for a comparison like this: How many 

times would I "have to run, say, Col or Myklestad to get the regression analysis as opposed to 
how many times I would have to run it if I Just hooked it up to an optimization program and Just 
reran it every time? Do you understand the question? 

Wei : I'll tell you. It depends upon how many design variables you are using. Right now we are 

using four independent variables. Normally we are using the quadratic regression analysis and 
here we have 36 cases. I personally believe that if we have less design variables and directly 
hook onto the analyzer combined with the optimizer, we are going to save time. If you have a 
tremendous [number of] design variables the regression analysis could be beneficial. I think 
the tradeoff here in independent variables is around seven; [this] would be a nice number. 

Bob Blackwell. Sikorsky Aircraft : I might ask if you could comment on whether the blade model 
and the inflow model and so forth that are used for your study are really sufficient for predic- 
tion of vibratory shears and prediction of blade response. Is it your [opinion] that a model 
as simple as this and able to be run for 36 times is sufficient or does the model have to get so 
detailed it just becomes cumbersome even with that? 

Wei : I personally feel that the present model still has to be improved so that we can use it 
for future design. Right now we only deal with four different independent variables and more 
independent variables are required in the future if we are going to do more in a real study. 
However, one thing that I can mention is that the people at Kaman [are] using the optimization 
technique to design for the SH-2 and they are using it now. How good are the results going to 
be? I don't have any answer at this moment. But we are going to see. 

Bob Taylor. Boeing Vertol : Just a quick question. Do you have any plans to do any testing to 

back up your theory? 

Wei : That's what I am saying. We are going to do the SH-2 composite rotor to hook on the SH-2 

helicopter. 

Taylor : That's how you are going to prove your theory? Build a full-scale blade? 

Wei : No, I can't give you an answer for that. 

Bob Goodman. Sikorsky Aircraft : It seems that the only way that you can really check this kind 

of thing is to run a variety of cases— isn't that true? I mean, really you need a baseline. 

Wei : We need ? data base to generate equations. I think, Bob, you can give more details. 

Bob Jones, Kaman : The regression equations are never going to be any better than the data 

base. If you have no faith in C8l then this is lousy. If you have no faith in something else 
then it is lousy. What you are doing is fitting statistical [variables to the] data base. If 
you have a good fit then it's a good equation, but it's no better than your data base, however. 
And you can do this with testing. I can get a data base with testing, fit a curve, [and do some 
interpolating]. This [fit] is really what it's based on. So there is no proof of theory if you 
want to look at it from that standpoint. We are worl^ g on methods where we have our regression 
equations based upon analysis and change them as we get testing results. 

Jing Ven, Bell Helicopter : John, I am Just curious to ask you what kind of inflow model you 

used here. 

Wei: We Just used the simple one that you see in the C81. 

Yen: You did not use the Dr. Gene Sadler's free wake [analysis]? 

Wei: No. 
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